DETERMINISTIC APPROXIMATIONS OF RANDOM REFLECTORS 

OMER ANGEL, KRZYSZTOF BURDZY AND SCOTT SHEFFIELD 

Abstract. Within classical optics, one may add microscopic "roughness" to a macroscop- 

ically flat mirror so that parallel rays of a given angle are reflected at different outgoing 

angles. Taking the limit (as the roughness becomes increasingly microscopic) one obtains 

^vj a flat surface that reflects randomly, i.e., the transition from incoming to outgoing ray is 

^-H described by a probability kernel (whose form depends on the nature of the microscopic 

^^ roughness). 

Cn We consider two-dimensional optics (a.k.a. billiards) and show that every random reflector 

5h on a line that satisfies a necessary measure-preservation condition (well established in the 

^T^ theory of billiards) can be approximated by deterministic reflectors in this way. 



I — 7 1. Introduction 

Ph This article addresses the question of which random reflectors can be approximated by 

piecewise smooth surfaces that reflect light according to the classical rule of specular re- 
flection which says that the angle of reflection is equal to the angle of incidence. One of 
'""' the basic results on billiards says that a certain measure on the space of pairs consisting of 

^ location and angle of reflection is preserved by every reflecting surface (see [T71 Thm. 3.1] 






or [3l Lemma 2.35]). Our main result, Theorem 2.3, shows that, except for this universal 



O 

00 restriction that applies to all reflecting surfaces, one can approximate in a weak sense every 

random reflector by a sequence of specularly (deterministically) reflecting surfaces. 

O A special case of this theorem is that one can approximate every deterministic reflector 

T-H that preserves the appropriate measure, including, for example, the reflector that reverses 

K>" the direction of each incoming ray. To see how counterintuitive this is, imagine a higher 

/S dimensional analog of this reflector: a mirror that appears entirely the color of the observer's 

H 
^ eyeball, because the only light rays traveling from the mirror to the eyeball are those that 

bounced off of the eyeball before reaching the mirror. This striking effect is implemented 

in practice with limited accuracy (for theoretical and practical reasons) in "retroreflectors" 

(see [13) and reflective paint (see [H]). 

There are several sources of inspiration for our project. The article [1] studies reflected 

Brownian motion with inert drift. In a follow up project, Z.-Q. Chen and the second author 

plan to study the limiting situation when the diffusion coefficient of reflected Brownian 
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motion goes to zero. It appears that in the hmit, the particle will move along straight lines 
with a random angle of reflection. 

The physics literature on random reflections is quite rich. We will not review it here; an 
excellent review can be found, for example, in [5]. 

Some early mathematical articles that considered random reflections were [m |9], and 
more recent ones include [H IS El [71 El III [I2l ttSl [16] . Many of these articles are concerned 
with the so called Knudsen reflection law. Knudsen proposed the cosine reflection law, in 
which the direction of the reflected molecule has a rotation-invariant distribution around 
the surface normal, with density proportional to the cosine of the angle with the normal, 
independent of the incidence direction. Knudsen's law is a model for gas dynamics. The 
same law is known as Lambert's cosine law in optics (see [2], pp. 147-148 or [13], Chap. 6). 

On the technical side, our results seem to be related, at least at the intuitive level, to the 
"digital sundial" theorem proved by Falconer ([ID], Thm. 6.9). Roughly speaking. Falconer's 
theorem says that there exists a set with prescribed projections in almost all directions. 

Although our article is close to the literature on billiards at the technical level, we will use 
the language of optics because our model is much closer to this circle of ideas at the intuitive 
level. 

We will state our results in a rigorous way in Section [2j Section [3] contains short proofs 
of the most elementary results. The proof of our main result. Theorem [2. 3[ is a multistage 
construction presented in Sections |4Jl5J 



2. The main result, conjectures and open problems 

For a very detailed and careful presentation of the billiards model in the plane see Chap. 2 
of [3]. We will be concerned with mirrors (walls of billiard tables) of very special shape. 
They are supposed to model macroscopically flat but rough reflecting surfaces. The paper 
of Feres [11] contains a rigorous mathematical presentation of this physical phenomenon and 
detailed analysis of its fundamental properties. Our setup is slightly different from that in 

Consider the following assumptions about a planar set M. These conditions contain, 
among other things. Assumptions A1-A4 from P^ Sections 2.1 and 2.4]. 

(Ml) M C {(xi,X2) eR'^:x2< 0}. 

(M2) For every k < oo, the set {(xi,X2) € M : —k < Xi < k} is the union of a finite 
number of compact C^ curves Tj. 

(M3) The curves Tj intersect only at their endpoints. Each curve either is a line segment 
or has non- vanishing curvature of one sign (it has no inflection points). The curves do not 
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form cusps at the intersection points, that is, the curves form an angle different from at 
the point where they meet (if there is such a point). 

We will say that M e Mi ii M satisfies (M1)-(M3). 

Let D = 'R X (— TT, 0) and define a cr-finite measure A on D hj A{dx, da) = —dx sin a da. 

A light ray can be represented as (x, a) = (x(t),a(t)), where x{t) is the location of the 
light ray at time t and a{t) G (— vr, vr] is the angle between the direction of the light ray 
and the positive horizontal half-axis, measured in the counterclockwise direction from the 
half-axis. Time will play no role in our arguments so it will be suppressed in the notation 
most of the time. We will always assume that light rays reflect from surfaces comprising 
M G A^i according to the rule of specular reflection, that is, the angle of incidence is equal 
to the angle of reflection, for every reflection. 

Let L* := {(xi,a;2) G M^ : X2 = 0}. It will be convenient to identify L^, with M, for 
example, we will consider A to be a measure on L^ x (— vr, 0). Consider the following natural 
condition. 

(M4) Suppose that M G M.i and for A-almost all (x, a) G -D, a light ray starting from 
(x, a) and reflecting from surfaces comprising M will return to L^, after a finite number of 
reflections. 

Condition (M4) is far from trivial; see, for example. Sec. 2.4 in [3] on accumulations of 
collisions. Moreover, some light rays reflecting from some mirror sets M E M.i will not 
return to L^,. We will show that (M4) holds for a large class of sets M . 

(M5) Suppose that M G M.i. Let {A'^}k>i be the family of all connected components of 
the open set {(xi,X2) G ]R^\M : X2 < 0}. There exists a subfamily {Aj^}k>i of {v4'^}fc>i (that 
is, every set Ak is equal to some set A^), such that every set Ak is bounded, L* C ljfc>i "^^fc, 
and the set {dAj n dA^ H L^,j,k > 1, j 7^ k} has no accumulation points in L^. 

We will say that M G A^2 if M G A^i and it satisfies (M4). 



ml3.1 Proposition 2.1. If M E Mi satisfies (M5) then it satisfies (M4) and, therefore, M E M 



2- 



Consider some M E M2- Suppose that a light ray starts from (xo,ao) with xq G L* 
and ao G (— 7r,0) at time 0, reflects from surfaces of M and returns to L^, at a time t, i.e., 
(x(t— ), a(t— )) = (xi,/3'), Xi G L^,, and t > is the smallest time with this property. Let 
/3 = /?' — vr. This defines a mapping K : D ^ D, given by K{x,a) = {y,(3). Clearly, K 
depends on M. 

We will write P(x, a; dy, d(5) to denote a Markov transition kernel on D, that is, for fixed 
(x, a) E D, P(x, a; dy, df3) is a probability measure on D. We assume that P satisfies the 
usual measurability conditions in all variables. 
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We will use Sx{y) to denote Dirac's delta function. Recall the transformation K and let 
Fk be defined by FK{x,a]dy,d(3) = SK{x,a){y, f3)dyd/3. In other words, Px represents a 
deterministic Markov kernel, with the atom at K{x,a). 

If /i„, n > 1, and /ioo are non- negative a-finite measures on some measurable space F then 
we will say that /i„ converge weakly to /ioo if there exists a sequence of sets Tj, j > 1, such 
that UiM Tj = r, fJ^niXj) < oo, fiociXj) < C)0 for all n and j, and for every fixed j, the 
sequence finiXj) converges weakly to fioo{Tj). 

j7. 1 Theorem 2.2. (i) Consider the transformation K : D ^ D corresponding to any M G A^2- 

The transformation K preserves measure A, that is, for any A G D with A{A) < oo, we 
have A{K^^{A)) = A{A). Moreover, K is "time reversible" in the sense that if K{Ai) = A2 
then K{A2) = Ai. 

(a) Suppose that for some sequence of sets Mn G A^2; corresponding transformations Kn, 
and some Markov transition kernels F{x,a;dy,dl3), we have 

A{dx,da)FK„{x,a;dy,d(3)^A{dx,da)F{x,a;dy,d(3) (2.1) ml6. 1 

in the sense of weak convergence on D^ as n -^ 00. Then P is symmetric with respect to A 
in the sense that for any smooth functions f and g on D with compact support we have 



f{y,l3)¥{x,a;dy,dl3)g{x,a)A{dx,da) = / g{y, /3)¥{x,a;dy,dl3)f{x,a)A{dx,da). 

(2.2) [3l2T" 
In particular, A is invariant in the sense that 

f{y,j3)F{x,a;dy,d/3)A{dx,da)= / f{x,a)A{dx,da). (2.3) ml2. 1 

D2 Jd 

See fn\ Sect. 4] for a similar result stated in a slightly different setting. The first part 
of the theorem says that all specular reflections are time reversible and preserve a certain 
measure. For this reason, A is known as the invariant measure for the collision map in the 
theory of billiards ([3], Sec. 2.12). This is related to Lambert's cosine law in optics (see [2], 
pp. 147-148 or [13], Chap. 6), also known as Knudsen's cosine reflection law in the context 
of gas dynamics (see |5j). The second part shows that this condition can be interpreted as 
symmetry for a Markov kernel (see (2.2)). This symmetry is preserved under weak limits of 
Markov kernels. The next theorem, which is our main result, says that the symmetry of the 



Markov kernel expressed in (2.2) is the only condition on a Markov kernel P necessary for 
the existence of deterministic approximations of random reflections represented by P. 

Recall that 6x{y) denotes Dirac's delta function. Suppose that the probability kernel P 
in Theorem 



2.2 



(ii) satisfies F{x,a;dy,df3) = 6x{y)dyF{x,a;df3) for some P. Heuristically, 
this means that the light ray released at x is instantaneously reflected from a mirror located 
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infinitesinially close to L^,. Then (|2.2|) and (|2.3|) imply that for all smooth bounded functions 

g{f3)F{x, a; d[5)f{a) sin a da, 



(-7r,0)2 



/ and g on (— tt, 0), and almost all x, 

/ f{l3)¥{x,a]d[5)g[a)smada 

i(-7r,0)2 

and 

/(/3)P(x, a; d(3) sin a da 

(-77,0)2 



(2.4) 



/(a) sin a da. 



(-TT.O) 



Theorem 2.3. Suppose that ¥{x,a;dy,d/3) = 6x{y)dy'P{x,a;d/3) where P satisfies (2.4). 
r/ien it/iere exists a sequence of sets Mn G M.2 o-nd corresponding transformations Kn such 
that 

A{dx, da)FK„{x, a; dy, d(3) — > A{dx, da)F{x, a; dy, d(3) 

weakly on D^ as n —t- oo. Moreover, Mn can be chosen in such a way that 

(a) Mn C {(xi,2;2) : —1/n < X2 < 0}, and 

(b) for every e > there exists D^ G D with A(D \ D^) < e such that all rays starting in 
Ds reflect from Mn exactly twice before returning to L* . 



Remark 2.4. The sets M„ that we construct in the proof of Theorem 2^ are not connected. 
We believe that the following conjecture can be proved using an iteration of the construction 



used in the proof of Theorem 2.3 



Conjecture 2.5. One can construct sets Mn so that they satisfy Theorem 2.3 except for (b), 
they are connected and every light ray reflects from Mn four times, except for a set of rays 
of A measure 1/n. 

The following problem is inspired by the "digital sundial" theorem of Falconer (see 
Thm. 6.9). 



Problem 2.6. Is it possible to construct sets Mn so that they satisfy Theorem 2^ and every 
light ray reflects from Mn only once, except for a set of rays of A measure 1/n? 



We believe that an analogue of Theorem Z3 holds in higher dimensions. 

3. Deterministic reflections 



Proof of Proposition 2. 1 . Assumptions A1-A4 from [3l Sections 2.1 and 2.4] are built into 
(M1)-(M3). Hence, there are no accumulation points for reflections from M, by the results 
in [31 Sect. 2.4]. 

The remaining part of the proof is based on [T71 Thm. 3.1 and Sect. 7.1]. We will only 
outline the main steps. Fix some k and consider Ak to be a bounded billiard table. We define 



ml2.2 



(2.5) ml2.3 



trans 
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the "billiard ball map" T as follows. Let v{x) be a continuous unit tangent vector field on 
dAk- For a light ray starting from a point x G dAk, let a be the angle between its direction 
and v{x). The light ray will hit dA^ at a point y G dA^. Let /3 be the angle formed by the 
direction of the light ray just after reflection at y and v{y). Then we let T{x,a) = {y,(3). 
Let dx represent the arc length measure on dAk. By p!7l Thm. 3.1], T preserves the measure 
dx sin ada on dA^ x (0,7r). Poincare's Recurrence Theorem (see [Tfl Thm. 7.4]) shows that 
{dx da)- almost all rays starting at {dAk fl L^:) x (0, vr) return to {dAk fl L*) x (0,7r) after a 
finite number of reflections. See pTf page 116] for more details. □ 



Proof of Theorem 2.2. (i) The claim that A{K ^{A)) = A{A) is a special case of a well 



known theorem, see [T71 Thm. 3.1] or [3] Lemma 2.35]. The fact that the light reflection 
process is time reversible implies that if K{Ai) = Ai then K{A'2) = Ai. 

(ii) Suppose that M„ G A^2 and Kn is the corresponding transformation. Consider some 
sets Ai,A2 G D and let f{x,a) = lAi{x,a) and g{x,a) = lA2{x,a). Then, using part (i), 

f{y, f3)FK^{x, a; dy, df3)g{x, a)A{dx, da) 
1ai {y, l3)5K„{x,a){dy, dl3)lA2 {x, a)A{dx, da) 

l^-i/^ Jx,a)lA2{x,a)A{dx,da) 
D2 " ^ '^ 

= A{K-\A,) n A,) = A{k;^\a,) n A,). 

For the same reason, 

g{y, /3)FkAx, a; dy, d/3)f{x, a)A{dx, da) = A{K-\A2) n Ai), 



lD2 

so 

f{y, P)Fk„{x, a; dy, d[5)g{x, a)A{dx, da) = / g{y, P)Fk„{x, a; dy, d[5)f{x, a)A{dx, da). 

D2 JZ)2 

A standard argument based on finite linear combinations of step functions and bounded 



convergence shows that (2.2) holds for smooth / and g with compact support and P of the 



form P_ft:„. This and weak convergence imply that (2.2) also holds for Markov kernels P that 



satisfy (2.1 ). We obtain (2.3) from (2.2) by the monotone convergence theorem applied to a 



sequence (7ri t 1 as n — > oo. □ 

4. Transposition reflector 

This section is devoted to the construction of a set of mirrors that transpose thin bundles 
of light of appropriate angular width. The first challenge is to guide bundles of light rays so 
that only an arbitrarily small amount of light is scattered in an unaccounted for way. The 
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J12.2 



second challenge (easier to present in a rigorous way than the first one) is to place various 
sets of mirrors so that they do not interfere with each other. 
Let A denote the closure of a set A. 



D ^ D a, simple symmetric function if there exists a 

X 



{x'l,x^) 



{a'l, 0^2) C D, k >1, such that 



Definition 4.1. We will call K 
countable family of rectangles Qk 
(i) Qk n Qj = for k ^ J, 
iii)[Jk>iQk = D, 

(iii) for any j, k, either {xi,X2) = {x{,X2) or (x^,X2) H {x{,X2) = k 
(iv) for every a < 00, there is only a finite number of k such that 
(v) for every k there exists a such that K{Qk) = {xi,X2) x {a}, 
(vi) for every k there exists j such that 

sin ada 



(x^,4)n(- 



-a, a 



^ 



sin ada, 



{x\,xX) 



{x\,x^^, K{Qk) C Qj and K{Qj) C Qk- 



We will call sup;.(x2 



•Ij 1 



V sup;,(a2 — a^) the mesh of K. 



Recall from Section [2] that for K : D ^ D, the kernel P;^ is defined by ¥k{x, a; dy, d(3) 



^K{x,a){y,P)dydl3. Suppose that P satisfies (2.4). Then standard arguments (see the proof of 
in Section™ show that there exists a sequence {Kn}n>i of simple symmetric 



Theorem 



2.3 



functions such that P^ (x, a; dy, d(3) — )■ 6x{y)dyF{x, a; d(3) weakly as n — )■ 00. Note that we 
do not claim that Kn arise as functions associated to reflecting sets in A^2- To prove Theorem 
it will suffice to show that for any simple symmetric function Kn : D ^^ D there exists 



2.3 



a sequence of sets M„ G M.21 n > 1, satisfying conditions (a) and (b) of Theorem 2.3 and 
such that if the functions Kn : D ^^ D correspond to M„'s then Kn — Kn -^ pointwise. 
The rest of the paper is devoted to the construction of sets M„. 

For p > 0, let Kp be the measure A restricted to Dp := {—p,p) x (— vr, 0) and note that 
the total mass of Ap is 4p. 

Lemma 4.2. For any Eq > there exist no < 00 and po > such that for any n > uq, 
p G (0, Po) o.iT'd any simple symmetric function K with mesh less than or equal to 1/n and 
such that {—p,p) is one of the intervals {x'l,X2), the following holds. 

(i) There exist a hounded set N G A^i and a set D* C Dp such that Ap{D*) < ApsQ. 

(a) The function Kn{x, a) corresponding to N is defined on Dp \ D*. 

(iii) For all (x, a) E Dp\ D^, we have \Kn{x, a) — K{x, a)\ < Eq. 

The function K^{x,a) is not necessarily defined for all (x, a) G Dp because some light 
rays starting from L^, and refiected in N may never come back to L*. 
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Proof of Lemma \4-^ We will omit some details of our arguments and estimates that are 
totally elementary but are tedious to write down. However, we will provide now a solid 
justification for the approximate formulas that will form the core of the proof. Suppose 
that two (parts of) ellipses are fixed and serve as mirrors for light rays. Suppose that a 
light ray starts from x & L^, at an angle a G (— 7r,0), then makes two reflections from 
the elliptic mirrors, and then returns to L^, at y G L^, and angle /3. An important (but 
elementary) observation is that the function {x, a) — )■ {y, (5) is analytic. The reason is that 
the equations that determine the points and angles of reflection of the light ray are quadratic 
with "parameters" that are analytic functions of x and a. A similar remark applies to other 
quantities that are functions of x and a, such as the location and angle of incidence when 
the light ray hits one of the ellipses, or the distance from the light ray to a fixed point. 

Step 1. In this step, we will construct a pair of elliptic mirrors that interchange two very 
thin bundles of light rays. This step is devoted only to the construction (or, in other words, 
definition) of the two mirrors. The next step will contain the proof that the mirrors reflect 
the light rays in the desired directions. 

We will use letters A, B and C (with subscripts and superscripts) to denote points in the 
plane. The notation AB will refer to a line segment with given endpoints; ABC will denote 
an arc of an ellipse passing through the three points. 

Consider some (small) ci > 0, a,/3 G (— vr + ci, — Ci) and (small) Aa, A/3 G (0,1), such 
that Aasina = A/3 sin /3. Suppose that p > and (see Fig. [l]), 

Ai = (-p,0), A2 = (0,0), A3 = (p,0). 

First assume that a = /3. In this case, let £a = £p be the arc of a circle represented in 
complex notation as £a = {re^^ : a — Aa < 7 < a + Aa} for some r > 0. Note that the 
circle is centered at A2. We will consider E^ to be a mirror. Light rays starting from (^,7) 
in D]^ := AiA^ x (a — AQ;,a; + Aa) will mostly return to the same set if p is small. More 
precisely, 

Vci > V771 > 3po < 00 Vp G (0,po) V« G (-TT + ci, -Ci) : 

Km "^'-'^ 

The case when a = (3 is rather easy so we will focus on the case a 7^ /3 in the rest of the 
proof. 



Recall notation from Definition 4.1 and let a^ = (aj + 02) f"^- Consider j and k such that 
K{Qk) C Qj and K{Qj) C Qk- If {a^ - a{) V {ai - a^) < 6o/2 then let K'{x, a) = (x, a'') 
for (x, a) G Qk and K'{x, a) = (x, a^) for (x, a) G Qj. Otherwise, we let K'{x, a) = K{x, a). 
Note that K' is a simple symmetric function and \K'{x, a) — K{x, a)\ < eo/2. Hence, it will 
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Figure 1. A single light ray interchange mirror system. 



figl 



suffice to find a set A^ and the corresponding function K]^ such that \Kisi{x, a) — K'{x, a)\ < 
£0/2. 

Consider the case when |a — /3| > £0/2 and a — Aa, a + Aa, (3 — A/3, (3 + A/3 G (— vr, 0). 
For any angle 7 G (— 7r,0), let L^ be the line passing through A2, with the slope tan(— 7). 
Unless stated otherwise, we will consider only points and sets below the line L*. 

Let B2 be a point on Lq,, and r^ := 1^2 — -82" |- Similarly, let -83 be a point on L/3, and 
r^ := |v42 — -B2 I (values of r^ and ry3 will be specified later). Let L' be the line passing 
through B2 and -B2 , and let C be the point on L' such that Si := \C — -Bf | = |C — i?2 |- 

Recall that for any ellipse, if a light ray leaves one of the foci and meets a point on that 
ellipse, it will refiect off the ellipse and pass through the other focus. 



Let 5" G La-Aa/2 and B^ G i^Q-+Aa/2 be points such that £a '■= B^B^B^ is an arc of an 
ellipse with the foci A2 and C (see Fig. l|. Similarly, let i?f G L ""'^ ^ 

be points such that Sp := -Bf i?2 -B3 is an arc of an ellipse with the foci A2 and C . 



/3-A/3/2 and i?3 G i^/3+A^/2 
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The rest of this step is devoted to informal estimates aimed at finding values of r^ and 
r^ that will define mirrors with desirable properties. Consider the bundle B of light rays 
starting at A2 along lines between La-Aa/2 and La+Aa/2- We will find parameters of our 
construction such that for the resulting mirrors £a and Sjs, "most" light rays in B reflect 
from £a, then reflect from £j3 and then hit A2. Let pi be the width of B close to £a, just 
before hitting Sa- Since Aa is assumed to be small, the rays in B are almost parallel and S^ 
is almost fiat. Hence, the width of B close to £a, just after hitting £a is pi(l + o(Aa)). We 
have pi = ra{Aa + o(Aa)) for small Aa. 

Since Sa is an arc of an ellipse with foci A2 and C, all light rays in B will pass through C 
and then they will hit Sf^, assuming that S13 is large enough. Since Sj^ is an arc of an ellipse 
with foci A2 and C, the light rays in B reflected from S13 will hit L* at A2. 

Let P2 be the width of B close to Sp, just before hitting S^. The width of B close to £^^, 
just after hitting S^ is p2{l + o{Aa)) . We want the bundle of light rays B to form a cone with 
angle A/3 at the point A2 on the way out. So we would like to have p2 = t^(A/3 + o(A/3)) 
for small A/3. The point C is half way between B2 and B2 so we would like to have pi = P2- 
Therefore, we choose Tq, and r/3 so that they satisfy 

r,(Aa + o(A«)) = p^ = p2 = r^(A/3 + o(A/3)). 

More precisely, we choose r^ and r^ so that 

- = ^ = ^- (4.1) [W:t 

rp Aa sin/3 1 

The above formula incorporates our previously made assumption that Aasina = A/3 sin /3. 

Step 2. We will now argue that a pair of mirrors defined in the previous step guides a thin 
bundle of light rays to the appropriate exit location and exit angle, with arbitrarily little 
"waste" , for appropriate values of parameters of the construction. Our argument is based on 
detailed analysis of small changes in the initial conditions (the starting location and angle 
of the light ray) on the "output," that is, the location and angle for the light ray exiting the 
lower half-plane. 

Let Ea be the ellipse of which £a is a part and let E^ have the analogous meaning. The 
informal estimates of Step 1 have the following rigorous version. For any ci, £ > there exists 
El G (0,e) such that for all a,/3 G (— vr + ci, — Ci) with |a — /3| > eo/2, Aa, A/3 G (0,£i), and 



Va and Tp satisfying (4.1), we have the following. If a light ray starts from ^2, follows La+a^ 
with «! G [—(1 — £i)A«/2, (1 — ei)/\a/2\ and reflects from £a then it will intersect E^ at a 
point X such that 



\x 



BP,\<{l-e){\B^,-Bl\y\B^,~B^,\). (4.2) \^21?2 



DETERMINISTIC APPROXIMATIONS OF RANDOM REFLECTORS 



11 



L 



A^ A^ A2 Bg A^ 




Figure 2. Trajectories of various light rays. 



fig2 



Similarly, if a light ray starts from A2, follows -^^+/3i with /3i G [— (1 —£1) A/3/2, (1 — ^i) A/3/2] 
and reflects from Ss then it will intersect E^ at a point y such that 



\y-B^\<il 



e){\B'^-B:^\y\B^-B!^\). 



(4.3) S21.3 



It follows from (4.2)-( [43 ) that all light rays starting from A2 at angles in [— (1— £i)Aa/2, (1 — 
£:i)Aa/2] or [—(1 — ei)A/3/2, (1 — ei)Af3/2] reflect from both mirrors Sa and Sp and then 
reach L^, again at A2. 

In the following part of the argument, we will consider a,P,ra and r/3 to be "fixed" 



parameters satisfying (4.1) and we will choose Aa, A/3 and p sufficiently small, so that 



certain conditions are satisfied. 

We will define a number of arcs, lines and points which will be used in the next part of 



the proof; see Fig. [2j Consider any point A4 G A1A3 and light ray -Ri emanating from A4 at 
an angle rj. Let -B4 G Eq, and Bq G E^ denote the points where this light ray refiects in E^ 
and E^. Let R4 be the line between B4 and Bq. Let Rq be the line that contains this light 
ray after refiection in E^ and let Bg be the point where Rq intersects L^,. Let R2 be the line 
passing through A2 and parallel to -Ri and let B^ be the intersection point of R2 with E^. 
Let R3 be the line passing through A2 and B4. Let R5 be the line that represents light ray 
i?3 after reflecting from Eq,. Let B-^ be the point where R5 intersects E^. Let R^ be the line 
passing through A2 and Bq. 



12 
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Let «! be the angle between A2B2 and the tangent to £a at -B^. The angle ai is a function 

of a, /3, Ta and r^ only. Assume that r/ G [a — (1 — ei)/S.a/2, a + (1 — £:i)Aa;/2]. Then 

\B, -B^\<{l-er + o(£i))(|i?r - B^\y \B^ - E^^j) 

r,(Aa/2)(l + o(Aa)) 



We also have 



and, therefore, 



< (l-£i+0(£i))- 



\B^-B^\<{l + o{p)) 



smcti 



psma 
sinai ' 



\B^-B^\<{l-e^ + o{e{)) 



r„(Aa/2)(l + o(Aa)) 



+ (l + o(p)) 



psma 



smai smai 

It follows that for any Ci,e > there exist ei G (0,£:), £2 > such that for all a, /3 G 
(— TT + ci, — ci) with |a — /3| > eo/2 and Aa, A/3 G (0,£:i), r^ and r/3 satisfying (4.1), all 

p G (0, e2(Aa A A/3)) and r^ G [« - (1 - £:i)A«/2, a + (1 - £:i)A«/2], 

we have -B4 G Sa- By symmetry, if r/ G [/3 — (1 — £i)A/3/2,/3 + (l — £:i)A/3/2] then Ri intersects 
E^ at a point in Sp. 

Suppose that rj E [a — (1 — £:i)Aa;/2, a + {l — £:i)Aa;/2]. Let /3i be the angle between A2B2 
and the tangent to £fs at Sg • Note that /3i is also the angle between between CB2 and the 
tangent to Sp at B2. Let 71 be the angle between Ri and R^. The angle between R4 and R5 
is also 7i. We have 

7i = sin a |y44 — ^21(1 + o{p) + o{Aa))/ra < psinQ;(l + o{p) + o{Aa))/ra- 

Thus 

2si sin7i 



\Br — Bj 



< 



sin/3i 
2sipsinQ;(l 



l + o(p) + o(A/3)) 
o(p) + o(A«)) 



(4.4) 



Tq sin /3i 



l + o(p) + o(A/3)) 



2sipsina 
Tq sin /3i 



o(p) + o(A/3)). 



It follows from (4.2), (4.3) and (4.4) that for any Ci,£: > there exist ei G {0,e), £2 > 
such that for all a, /3 G (— tt + ci, — ci) with |a — /3| > £0/2 and Aa, A/3 G (0, ei), Tq, and r/3 
satisfying (4.1 ), and all p G (0, £:2(Aa A A/3)) we have the following. If a light ray starts from 
a point in AiA^ at an angle 77 G [a — (1 — ei)Aa/2, a + (1 — £:i)A«/2] then it reflects from 
Sa and then it intersects S13. Similarly, if a light ray starts from a point in AiA-^ at an angle 
r] G [/3 — (1 — £:i)A/3/2, /3 + (1 — £:i)A/3/2] then it reflects from Sp and then it intersects £a- 
Let a denote the distance between C and R4. We have 



S21.1 



Sisin7i(l + o{p)) = sisina IA4 - ^21(1 + o{p) + o{Aa))/ra- 
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Note that if a light ray moved along CBq and reflected from E^ then it stays on R^. The 
angle between Rq and Rj is equal to the angle between CBq and R4 and, therefore, it is 
equal to 72 := (a/si)(l + o{p) + o{Aa)). The following estimate for the distance between A2 



and Bg uses (4.1) 



1^8 - A2I = ^(1 + o(A«)) = ^^(1 + o{p) + o(A«)) 
smp sismp 

Sisina|A4- Aalr/j 

= ^-5 (1 + o(p) + o(Att)) 

raSi sm p 

= |A4-v42|(l + o(p) + o(Aa)). 

We combine this estimate with the conclusions obtained so far to see that for any ci,e > 
there exist ei G (0, e), £2 > such that for all a, /3 G (— vr + ci, — Ci) with |q; — /3| > eo/2 and 



Aa, A/3 G (0,£i), r^ and r^? satisfying ( [4l| , all p G (0,e2(Aa A A/3)), and all A4 G A1A3 
satisfying |y44 — v42| < (1 — £2)^ we have the following. If a light ray starts from A^ at an 
angle rj E [a — (1 — £1) Aa/2, a + (1 — £1) Aa/2] then it reflects from Sa, next it reflects from 
Si3 and then it intersects L* between Ai and ^3. Similarly, if a light ray starts from A4 at 
an angle 77 G [/3 — (1 — ei) A/3/2, /3 + (1 — ei) A/3/2] then it reflects from S/3, next it reflects 
from £a and then it intersects L* between Ai and ^3. 

Step 3. We will now assemble a finite family of pairs of mirrors so that they properly guide 
light rays entering the system at "most" angles and the mirrors do not interfere with one 
another. The lack of interference will be achieved by inductive scaling of the mirrors, that 
is, making them large, so that a light ray traveling between a pair of mirrors takes a path 
far beyond all the mirrors constructed earlier in the inductive procedure. 

Recall the notation from the statement of the lemma. Let ci > be so small that 

Ap {Dp \ ((-p, p) X (-vr + ci, -ci)) ) < Apso/IQ. 

Note that Ci satisfying this condition can be chosen independently of p. 

Recall the simple symmetric function K' defined in Step 1 and let {Q'k}k>i = {(3^1,3^2) ^ 



(Q;^,a2)}fc>i be the family of rectangles as in Definition 4.1 Let J be the family of pairs 
(j, k) such that j < k, 

K'{{-p,p) X {a{,ai)) C {-p,p) x K,a^) 

and, therefore, 

K'{{-p,p) X (at, "2)) C {-p,p) X {ai,ai). 
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Note that according to Definition 4.1 (iv), the set J is finite. Let J be the set of all j such 
that \a\, a^] n (-vr + Ci, -Ci) 7^ 0. For (j, k) G J, let 

oP = {a{ + ai)/2, 
AS"' = al — a{, 
Aa'' = AS-' sin a-'/ sin a'^. 

Note that, typically, it is not true that AS'^ = 0:2—0;^ because if (j, fc) G J and j ^ k then 
(fc, j) ^ J. The quantity AS'^ is well defined because if (j, /c) G J and j ^ k then there is 
no i such that {k,i) G J. It is easy to see that for fixed Ci,ei > we can choose uq so large 
that for all n > no, if the mesh of K is smaller than 1/n then for all j E J we have 

[& - (1 - £i)ASV2, a^ + (1 - ei)ASV2] c («{, a^s)- 

Let 

e(ei) = (-TT + ci, -ci) n IJ p - (1 - ei)AaV2, S^' + (1 - ei)AaV2]. 

For a fixed ci > 0, we can make ei > smaller and uq larger, if necessary, so that 

A,((-p,p)xe(£i))>4p(l-£o/8). 

We make tiq larger, if necessary, so that for n > no and all j E J' we have AS-' < Si. 
Choose £3 > so small that 

Ap((-(1 - es)p, (1 - e;)p) X &{ei)) > 4p(l - bq/A). 

We will write r{a) instead of r^ and S{a) instead of Sa, for typographical reasons. 

Let A*a = miuj-gj^ AS-'. It follows from what we have shown in Step 2 that for ci,£: > 
there exist no, ^i G {0,e), £2 > such that for all p G {0,e2A*a), n > uq and all j E J' we 
have the following. If (j. A;) G J and a light ray starts from a point in {—{l — es)p, (1 — e3)p) C 
L^, at an angle 

r] G [S^' - (1 - £i)ASV2, a^ + (1 - £i)AS-'72] 

then it reflects from S{a^), next it reflects from S{a^) and then it intersects L* at a point in 
(— (1— e3)p, (1— £:3)p). Moreover, if a light ray starts from a point in (— (1 — £:3)p, (1 — £:3)p) C 
L^ at an angle 

r] G [S'= - (1 - £i)ASV2, S'^ + (1 - £i)ASV2] 
then it reflects from £{a^), next it reflects from S{a^) and then it intersects L* at a point in 

(-(1-£3)P,(1-^3)P). 

Let (ji, /ci), (j2, ^2), • • • be an arbitrary ordering of pairs in J. Choose r(S-'^) and r{a''^) so 
that they satisfy (4.1 ). Suppose that rla^"') and r{a'''") have been chosen for m = 1,2, . . . ,i. 
We choose r(S-''+^) and r{a''^+'^) so that they satisfy (|4.1) and they are so large that lines 
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connecting any points in S{&'+^) and S{a'^'+'^) do not intersect any £^(a-''") and £^(a'^™) for 
m = 1,2, . . . ,i. 

Recall that J is finite and let ni be its cardinality. For (j, k) G J, let Dp^y^^-^ be the set of 
light rays in Dp such that the ray starts from a point in [—(1 — e3)p, (1 — £3)p] C L* at an 
angle 

r] G [& - (1 - ei)ASV2, a^ + (1 - £i)Aa^/2] 



or at an angle 



77 e [S^ - (1 - £i)ASV2,a^ + (1 - £i)AaV2]. 



Let Dp^i^j^k) be the set of light rays in Dp^(^j^k) such that at some time the ray reflects from a 
set £{a"^), for some m 7^ j, k. li p = 0, that is, if [— (1 — e3)p, (1 — £:3)p] is the single point A2, 
then Dp^i^j^k) = 0, by the claim made in the previous paragraph. The proportion (in terms of 
the measure Ap) of light rays in Z^p (j_fc) which belong to Dp^(j^fc) is a continuous function of 
p > 0, with zero limit when p 10, that is, limp^^o Ap(Dp,(j,A:))/Ap(-Dp,(j,fc)) = 0. It follows that 
for sufficiently small p > 0, we have Ap{Dp(^jk)) < 4peo/(16ni). Let Dp = IJ(ifc)eJ -^p.O'.'^O- 
We see that for fixed r{a^)'s, we can make p > so small that Ap{Dp) < 4peo/16. 

Let A^ = [Jj^j^{<^^)- Our arguments have shown that the function Kj^ satisfies the 
assertions listed in the lemma, with function K' in place of K, and eo/2 in place of Eq. We 
have observed in Step 1 that this is sufficient for the proof of the lemma. □ 



sec: cantor 



5. Light reflector cells in Cantor set holes 

The system of mirrors constructed in the previous section is very "inefficient" in that it 
has a large diameter compared to the length of the segment {—p,p) of L* where the light 
enters and exits the system. In this section, we will build a "compact" version of the reflector 
using scaling and a Cantor-like construction. The following result is very similar to Lemma 



4.2 except that we use smaller mirrors. 



J12.3 



Lemma 5.1. Fix any p, pi > 0, a G M and let Ap be the measure A restricted to Dp : = 
{a — p,a + p) X (— 7r,0). The total mass of Ap is Ap. For any ^o > there exists tiq < 00 
such that for any n > uq and any simple symmetric function K with parameter n such that 
[a — p,a + p] is one of the intervals [x\, X2], the following holds. 

(i) There exist a compact set N E Aii and a set D* C Dp such that Ap{D*) < 4peo- 

(a) The function K]\f{x, a) corresponding to N is defined on Dp \ D* . 

(Hi) For all (x, a) E Dp\ D*p, we have \Kff{x, a) — K{x, a)\ < Eq. 

(iv) N C {a - p, a + p) X (-pi,0). 
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Figure 3. Sets M-l^ of three "generations." The sets M^ are convex by defi- 
nition. They are pictured as circles but in fact they are not discs. f ig7 



Proof. By translation invariance, it is enough to discuss the case a = 0. 

Recall that we consider {—p,p) to be a subset of the horizontal axis L^,. Let A^o be a 



new name for the set satisfying Lemma 4.2 and let M be the closure of the convex hull of 



A'o U (— p, p). Since Ci in the proof of Lemma 4.2 is strictly positive, M fl L^, = [— p, p]. Let 



ro be the diameter of M and let Mr = {x eM."^ : 3y E M such that x = {r/ro)y}, i.e., Mr is 
a dilation of M. 

We will define sets M^. Each of these sets will be a horizontal shift of Mr, that is, 
M/ = Mr+{b, 0) for some 6 G M depending on j and r. For any such set, we let M/ = M^nL^:. 
Note that the length of M^ is 2(r/ro)p. 

The sets M^ will be grouped into a countable number of families, with all sets in one family 
having the same size. We will pack smaller sets among the bigger sets in a tight way, as 
much as possible. The intersections of the sets M^ with L* will form a pattern qualitatively 
similar to holes in the Cantor set. See Fig. [3j 

For r G (0, 1), let Ur be the maximum number k of disjoint sets Mj^M'^, . . . , M^ such that 
Ui<j<fc ^r '^ {(^i5 ^2) : \xi\ < 1}. It is easy to see that for some C2, r^ > and all r G (0, r^,) 
we have rnr > C2- It follows that for some C3 > 0, for every r G (0, r*) one can find a family 
of disjoint sets M^,M^, . . . ,M^ such that IJi<i<fc ^r ^ {(a;i)3;2) : la^il < 1} and the total 
length of IJi<i<fc ^r is greater than C3. 

Choose ri G (0,pi) and a family of disjoint sets {^ri}i<i<ii such that IJi<j<ji ^n <^ 
{(a;i,a;2) : |xi| < p} and the total length of IJi<7<ji ^n is greater than c^p. 

We proceed by induction. Let Ai be the Lebesgue measure of the set [— p, p] \ IJi<7<ii ^n- 
We find r2 G (0,ri) and a family of sets {M^^}i<j<j2 such that 

(i) Ui<,<,2^r'2 C {(a;i,a;2) : ki| < p}, 
(ii) all sets in the family {M^^} k=i,2;i< j<jk are disjoint, and 
(iii) the total length of IJi<i<i2 ^4 is greater than C3A1/2. 
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The general inductive step is the following. Suppose that we have defined r^ > and 
families {M}^}i<j<j^ ior k = 1,2, . . . ,i. Let Aj be the Lebesgue measure of the set [— p, p] \ 
Ui<A:<i Ui<j<jfe ^rfc- We find Tj+i G (0,ri) and a family of sets {M}^^^}i<j<j^^^ such that 

W Ui<j<j,+i ^i^,+i C {(xi,X2) : |xi| < p}, 
(ii) all sets in the family {M^^}i<k<i+i;i<j<jk are disjoint, and 
(iii) the total length of IJi<j<ii i ^r, i is greater than c^\i/2. 

Let T/ be the linear transformation that maps M onto M^ and let A^^ = T^{Nq). Let 
^* ~ Ufc>i UkjXj, -^4' ■'■^ i^ ^^^y ^° ^^^ ^^^^ ^^^ Lebesgue measure of the set [—p,p] \ 
[Jk>i Ui<j<jfc ^4 i^ ^^^°- ■'■^ i^ ^^^° clear that A^* C {(xi,X2) : |xi| < p}. 

Let -D C -Dp be such that Ap{ID) > 4p(l — ^o) and for all (x, a) G -D, we have \K]\ff^{x, a) — 
K{x, a)\ < eq. Let Br = {{x, a) : 3{y, a) e 3 such that x = {r/ro)y}. If M,} = Mr + (6^, 0) 
then we let IDi = {{x, a) : 3{y, a) & D such that x = y + {pl., 0)}. Let A;^ be the measure A 
restricted to {—{r/rQ)p + b^, {r /rQ)p + 6^) x (— vr, 0). By the scaling and shift invariance of 
reflections, we see that for every r > and j we have Ai{Dj.) > (r/ro)4p(l — Eq) and for all 
{x,a) G D^, we have \KN^{x,a) - K{x,a)\ < Eq. Let D* = Df,\ |J^>-^ Uikjkj^ ^j.^- Then 
Ap{D*) < 4:pEo, and for all (x, a) G -Dp \ -D*, we have \Kn^{x, a) — K{x, a)\ < Eq. 

The set A^* satisfies all assertions in the lemma except that it does not belong to A^i 
because it consists of an infinite number of curves (arcs of ellipses) and hence it does not 
satisfy condition (M2) of Sectionpj To address this problem, we let A^ = IJi<fc<fco Ui<j<i ^4 
for some k^ < oo. It is easy to see that A^ and the accompanying function K^ satisfy all the 
conditions stated in the lemma if k^ is sufficiently large. □ 



Proof of Theorem 2.3[ Recall the definition of a simple symmetric function from the begin- 
ning of Section 4 Suppose that P satisfies (2.4). Recall that j_ sin ada = —2, let m > 1 



be an integer and let 7^ be defined by — JJ'' sin ada = fc2 "* for k = 0,1, . . . , 2"^+^. If we 



take /(a) = l(^m^^™ ](a) and g{a) = l(-^™^^m ](a) then (2.4) yields 



/ 



/ P(x, a; (i/3) sin arfa = / / P(x, a; (i/3) sin arfa. (5.1) o2.1 



' (7f ,7™n] J {7^,7,'r+i] -^ (7r>7r+i] -^ (7f ,7j;i] 

For a fixed integer — 00 < /c < 00 let 



a{n,i) = — I I j P(x, a; (i/3) sin adadx. 

'[fc/m,(fc+l)/m] i(7r.7r-n] >^(7f ,7™,.i] 



It follows from (5.1) that a{n,j) = a{j,n). Let /3™o = 7^ and let /3™j be defined by 

'(fc/m,(fc+l)/m] J7- 0<i<i-l 



/ / sin adadx = >, 

J(k/m,(k+l)/m\ J^:^ n^„-^» 
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for z = 1, . . . , 2*"+^ Note that /^_^ ^j F{x, a; dP) = 1 so 

2_] ^(n,^) = — / / sin adadx, 

0<i<2-+i-l J[k/m,{k+l)/m] ^(7^,7^+1! 

and, therefore, /3™2™+i = 7^+1- 

Let ^™ = (/3- + /3™+i)/2 and ir^(x,a) = (x,^- ) for (a;,a) G ik/m,{k + l)/m] x 
(/3j>,/3j>+i]- Since a(n, j) = a(j,n), we have 

sin ada = I sin ada. 

/om om ) 



4.1 



This shows that Km{x,a) is a simple symmetric function according to Definition 

Recall that ioi K : D ^i- D, the kernel Fk is defined by fxix, a; dy, df3) = SK(x,a){y, (3)dyd(3. 
It is elementary to check that F^ {x, a; dy, df3) — )■ 6x{y)dyF{x, a; df3) weakly on D as n — t- oo. 
Note that we do not claim that Kn arise as functions associated to reflecting sets M G A^2- 



To prove Theorem 2.3 , it will suffice to show that there exists a sequence of sets M„ G A^25 



n > 1, satisfying conditions (a) and (b) of Theorem 2.3 and such that if the functions 



Kn : D ^ D correspond to M„'s then K^ — K^ — )■ pointwise. 

Fix some Kn and let {x\,X2) be one of the corresponding intervals as in Definition 
Let Dk = {x\,X2) X (— 7r,0). According to Lemma 5.1 and its proof, there exist a compact 



4.1 



set Nk G A^i and a set Dl C Dk such that A{Dl) < eo2 '', and the following assertions hold. 

(i) The function K^^{x, a) corresponding to A^^ is defined on Dj^ \ Dl. 
(ii) For all (x, a) & Dk\ Dl, we have \KNi_{x,a) — Kn{x,a)\ < Eq, and the light ray with 

the starting position and direction represented by {x, a) refiects exactly twice from A^^ 

before returning to L*. 
(iii) NkC{xlxl)x{-l/n,0). 

Note that for j ^ k, either the sets Nj and A^^ do not intersect or they are identical. 



Moreover, the construction presented in Lemma |5.1| allows us to assume that light rays that 
start in D^ \ Dl do not leave (x^,X2) x (— oo, 0). In this sense, the families of mirrors A^^ do 
not interfere with one another. Let M* = IJfc>i ^^ ^^^ ^^^ ^ni^^^^) be the function on D 
corresponding to M*. Let D* = IJfc>i-^fc ^^'^ ^'^^^ ^^at A{D*) < Eq. We have 

(i) The function Kn{x, a) is defined on D\D*. 
(ii) For all {x,a) E D\ D* , we have \Kn{x,a) — Kn{x,a)\ < Eq, and the light ray with 

the starting position and direction represented by (x, a) reflects exactly twice from M* 

before returning to L^. 
(iii) M* cMx (-l/n,0). 
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It is easy to see that M* G A^i for every n but it is not necessarily true (actually unlikely) 
that M* E A4.2- We can instead consider sets 

M„ = M* U {(xi, X2) : X2 = -l/n} U |J{(xi, X2) : Xi = x'l or x^, -l/n < X2 < 0}. 

k 

It is easy to see that functions Kn corresponding to M„ satisfy all properties (i)-(iii) listed 
above. Moreover, M„ satisfy condition (M5) of Section [2] so M„ G A^25 according to Propo- 



sition 2.1 This completes the proof of the theorem. □ 



Acknowledgment. We are grateful to Persi Diaconis and John Sylvester for the most 
helpful advice. We thank the anonymous referee for many suggestions for improvement. 

References 

[1] K. Burdzy, R. Bass, Z. Chen and M. Hairer, Stationary distributions for diffusions with inert drift 

Probab. Theory Rel. Fields 146 (2010) 1-47. 
[2] S. Chandrasekhar, Radiative transfer. Dover PubUcations, Inc., New York, 1960. 
[3] N. Chernov and R. Markarian, Chaotic billiards. Mathematical Surveys and Monographs, 127. Amer- 
ican Mathematical Society, Providence, RI, 2006. 
[4] F. Comets and S. Popov, Ballistic regime for random walks in random environment with unbounded 

jumps and Knudsen billiards Ann. Inst. H. Poincare Probab. Statist. (2012) (to appear). 
[5] F. Comets, S. Popov, G. Schiitz and M. Vachkovskaia, Billiards in a general domain with random 

reflections Arch. Rat. Mech. Anal. 191, (2009) 497-537. 
[6] F. Comets, S. Popov, G. Schiitz and M. Vachkovskaia, Erratmn: Billiards in a general domain with 

random reflections Arch. Ration. Mech. Anal. 193, (2009), 737-738. 
[7] F. Comets, S. Popov, G. Schiitz and M. Vachkovskaia, Quenched invariance principle for the Knudsen 

stochastic billiard in a random tube Ann. Probab. 38 (2010) 1019-1061. 
[8] F. Comets, S. Popov, G. Schiitz and M. Vachkovskaia, Knudsen gas in a finite random tube: transport 

diffusion and first passage properties. J. Stat. Phys. 140 (2010) 948-984. 
[9] S. Evans, Stochastic bilhards on general tables. Ann. Appl. Probab. 11, 419-437 (2001). 
[10] K. Falconer, Fractal geometry. Mathematical foundations and applications. Second edition. John Wiley 

& Sons, Inc., Hoboken, NJ, 1990. 
[11] R. Feres, Random walks derived from billiards. Dynamics, ergodic theory, and geometry, 179-222, 

Math. Sci. Res. Inst. Publ., 54, Cambridge Univ. Press, Cambridge, 2007. 
[12] R. Feres and H.-K. Zhang, The spectrum of the billiard Laplacian of a family of random billiards. J. 

Stat. Phys. 141, (2010) 1039-1054. 
[13] R. Klette, K. Schliins and A. Koschan, Computer vision. Three-dimensional data from images. 

Springer- Verlag Singapore, Singapore, 1998. 
[14] S. Lalley and H. Robbins, Stochastic search in a convex region. Probab. Theory Relat. Fields 77, 

99-116 (1988). 
[15] M. Lapidus and R. Niemeyer, Towards the Koch snowflake fractal billiard: computer experiments 
and mathematical conjectures. Cems in experimental mathematics, 231-263, Contemp. Math., 517, 
Amer. Math. Soc, Providence, RI, 2010. 
[16] M. Lapidus and R. Niemeyer, Families of Periodic Orbits of the Koch Snowfiake Fractal Billiard. Math 

ArXiv 1105.0737 
[17] S. Tabachnikov, Geometry and billiards. Student Mathematical Library, 30. American Mathematical 
Society, Providence, RI; Mathematics Advanced Study Semesters, University Park, PA, 2005. 



20 OMER ANGEL, KRZYSZTOF BURDZY AND SCOTT SHEFFIELD 

[18] Wikipedia, Retroreflector, jhttp : //en . wikipedia . org/wiki/Retroref lector Online; accessed 4- 

March-2012. 

[19] Wikipedia, Scotchlite, http : //en. wikipedia. org/wiki/Scotchlite Online; accessed 4-March-2012. 

Department of Mathematics, University of British Columbia, 121 - 1984 Mathematics Rd, 
Vancouver, BC, V6T 1Z2, Canada 
E-mail address: angelOmath . ubc . ca 

Department of Mathematics, Box 354350, University of Washington, Seattle, WA 98195, 
USA 

E-mail address: burdzy@math.washington.edu 

Department of Mathematics, Massachusetts Institute of Technology, 2-180, 77 Mas- 
sachusetts Ave., Cambridge, MA 02139, USA 
E-mail address: sheffield@math.mit.edu 



